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Abstract 

We compute the next-to-leading logarithmic (NLL) ultrasoft running of the spin-independent sin- 
glet potentials up to 0(l/m 2 ), and the corresponding contribution to the spectrum. This includes 
the static energy at next-to-next-to-next-to-leading logarithmic (NNNLL) order. As a byprod- 
uct of these results we set the stage for the complete analytic and numerical computation of the 
heavy quarkonium spectrum with N 3 LL accuracy for I ^ (angular momentum) and s = (spin) 
states. We also compute the next-to-next-to-next-to-next-to-leading order (N 4 LO) ultrasoft spin- 
independent contribution to the heavy quarkonium mass and static energy. 

PACS numbers: 12.38. Cy, 12.38.Bx, 12.39.Hg, 11. 10. St 
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I. INTRODUCTION 



0,6 



The evaluation of the three-loop soft contribution to the static potential |2j has given 
the final piece needed for the complete evaluation of the heavy quarkonium spectrum at 
NNNLO 3, J]. This result is also necessary for the long term project of obtaining the heavy 
quarkonium spectrum with NNNLL accuracy 1 , and the non-relativistic sum rules and t-i 
production near threshold with NNLL /NNNLO one. It is then timely to compute the next- 
to-leading ultrasoft running of the potentials, as it enters in the evaluation of the previously 
mentioned observables with such precision. Therefore, in this paper we compute the next- 
to-leading ultrasoft running of the spin-independent potentials. For the static potential the 
running has already been computed in Ref. 7|, we confirm this result, for the 1/m and 1/m 2 
potentials the result is new. 

We will write the results for the potentials in position space. This could be eventually 
convenient for future numerical evaluations of decays and sum rules. Moreover, such compu- 
tation is interesting on its own, as it shows the subtleties appearing in the explicit matching 
(and the associated scheme dependence) between the soft and ultrasoft computation. In 
momentum space this discussion has already been made in Ref. 3j], we do so in position 
space (see also Ref. js] for a similar discussion for the QCD static potential in three dimen- 
sions, and Ref. {9] for the SUSY QCD static potential). This discussion is relevant, as the 
subtraction scheme of the ultrasoft divergences could be different in position and momentum 
space. 

The basis of potentials that we use for the 1/m and 1/m 2 potentials is redundant. There- 
fore, the expression obtained for each potential is ambiguous, since field redefinitions can 
shift some contributions from the 1/m to the 1/m 2 potentials, and viceversa. In principle, 
one could avoid this problem by considering all the 1/m and 1/m 2 potentials as a whole. Yet, 
this observation has implications in one of the possible applications of our result: compari- 
son (at short distances) with the recent lattice simulations of the 1/m and 1/m 2 potentials 



obtained in Refs. 10|, QJ| using their non-perturbative expression in terms of Wilson loops 
I2I . 13]. Unfortunately, the ambiguity just mentioned makes not possible a direct connection 
between the non-perturbative expressions of the potentials in terms of Wilson loops and the 



1 Actually, this accuracy has already been achieved for the hypernne splitting 0, Q . 
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perturbative computation. This would require a dedicated study that goes beyond the aim 
of this work. 

Whereas there are still some pieces left for a complete NNNLL evaluation of the heavy 
quarkonium spectrum, our result provides the missing link for the complete result for I ^ 
and s = states, the structure of which is given in this paper for the first time. The full 
explicit analytic form will be presented elsewhere. Another by-product of our computation 
is the N 4 LO ultrasoft spin-independent contribution to the heavy quarkonium mass. 

The next-to-leading ultrasoft running of the 1/m 2 potential was computed in Ref. [l^] 



in a different framework named vNRQCD (see Refs. |15l-ll7j). The computation of such 
object alone does not make much sense due to field redefinitions ambiguities that can shift 
contributions among different potentials. If we compare expressions for the 1/m 2 potential, 

(2) 

we disagree with the running of V^ 2 obtained in that paper but yet, as mentioned, this 
could be due to field redefinition ambiguities. In Ref. [18| the 1/m 2 and 1/m potentials 
were considered. Nevertheless, no definite outcome for the complete NLL ultrasoft running 
of the potentials was obtained. The reason was that the result was dependent on how the 
infrared divergences were regulated, so they even ended up having more than one possible 
result. Thus, a complete comparison with their results at this stage is not possible. Other 
issues that complicate the comparison are that the soft contribution is not included in their 
potential and the intrinsic scheme dependence. Either way, we believe that such future 
comparison should be better performed for specific observables. In this respect a future 
evaluation of the heavy quarkonium spectrum with NNNLL accuracy for I ^ and s = 
states with vNRQCD would be a good object for such comparison. 

Non-perturbative effects will not be considered in this paper and ultrasoft effects will be 
computed within perturbation theory (i.e. in the ma 2 ^> Aqcd limit), yet the renormal- 
ization group (RG) results will also be valid when ma 2 ~ Aqcd- The ultrasoft effects can 
be easily obtained in the non-equal mass case, as they only depend on the reduced mass. 
On the other hand the full soft contribution is only known in the equal mass case at the 
appropriate order. 

The outline of the paper is as follows. In Sec. [TT] we introduce the theoretical setup and 
show some bare results relevant for our computation. In Sec. [1111 we compute the ultrasoft 
correction to the spin-independent Hamiltonian with NLL precision and give explicit expres- 
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sions for the relevant potentials. In Sec. [IV] we give expressions for the heavy quarkonium 
energy and static potential with NNNLL precision (for the heavy quarkonium mass only 
when I 7^ and s — 0), and in Sec. |V]we present our conclusions. Finally, in the Appendix 
we gather some constants that appear throughout the computation. 



II. PNRQCD 

Up to NLO in the multipole expansion the effective Lagrangian density of pNRQCD takes 
the form QEJ: 



C us = Tr|s+ {id - h s {r)) S + 0* {iD - h {r)) O j 



+#V A (r)Tr{O t r-ES + S t r-EO} + g^^Ti {0 ] {r ■ E,0}} . (1) 

We define color singlet and octet fields for the quark- ant iquark system by S = S(r, R, t) and 
Qa _ Qa^ r ^ ^ respectively. R = m ™£ m x i + mi ™ 2 m2 x 2 is the center position of the system, 
and r = xi — X2. In order for S and O a to have the proper free- field normalization in color 
space they are related to the fields in Eq. (CQ) as follows: 

1 T a 

r -.S, = -^=O a . (2) 



All gluon and scalar fields in Eq. (CQ) are evaluated in R and the time t, in particular the 
chromoelectric field E = E(R, t) and the ultrasoft covariant derivative iD^O = idoO — 
g[A (R,t),O}. 

h s can be split in the kinetic term and the potential: 

P 2 

h s (r, p, Si, S 2 ) = h V s (r, p, Si, S 2 ), (3) 

2 m r 

h a (r, p, Si, S 2 ) = + V (t, p, S 1; S 2 ), (4) 

where m r = mim 2 / '(mi + m 2 ), p = — zV r and Si and S 2 are the spin of the quark and 
the antiquark respectively. For the equal mass case: mi = m 2 = m, the potential has the 
following structure (we drop the labels s and o for the singlet and octet, which have to be 
understood): 



,+•••, (5) 

m m A 
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v® = v® + vg\ 

= ^{P^^ 2) W} + ^L 2 + K (2) (r), (6) 
= V%j(r)L • S + V«(r)S 2 + ^(r)S 12 (f), (7) 

where S = Si + S2, Si2(f ) = 4(3r • Si r ■ S2 — Si ■ S2), and L = r x p. The last two equalities 
hold in 4 dimensions. L 2 is generalized to (q = p — p') 

iTxr 6 \ q J 

to be compatible with D dimensional calculations in momentum space. In this paper we 
focus on the ultrasoft corrections to the spin-independent singlet potentials. Therefore, we 



will not consider Eq. (J7J) (nor its generalization to D dimensions, see for instance Ref. 2l| ) 
in the following and drop the labels s and o, except for the static potentials. 

From now on we will use the index "5" to explicitly denote bare quantities. Parameters 
without this index are understood to be renormalized. The bare parameters of the theory 
are as [qb) and the potentials, generically denoted by Vb- o> and V{ Si0t A,B}(r) are the Wilson 
coefficients of the effective Lagrangian. They are fixed at a scale v smaller than (or similar to) 
1/r and larger than the ultrasoft and any other scale in the problem by matching pNRQCD 
and NRQCD. 

In our convention as is dimensionless and related to gs by (D = d + 1 = 4 + 2e) 

where v is the renormalization scale. It has a special status since it does not receive correc- 
tions from other Wilson coefficients of the effective theory. We renormalize it multiplicatively: 



where 



The RG equation of a is 



«b = Z a a , (10) 



00 1 

Za = l+5>«-. (11) 



-1 6 



v^—ot = aj3(a; e) = 2ta + a/3(a; 0) . (12) 
dv 
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In the limit e — > 

d d 
v—a = a/3(a; 0) ee a/3(a) = -2a— , (13) 
az/ da 

where 

^ = ^A. + • • • aft*) = -2a + + ••■), (14) 

and expressions for (3 , (3\, etc, can be found in the Appendix. 

The bare potentials Vb in position space have integer mass dimensions (note that this is 
not so in momentum space) and, due to the structure of the theory, we do not renormalize 



them multiplicatively (see the discussion in Ref. |22j). We define 



V B = V + 5V. (15) 

5V will generally depend on the (matching) coefficients of the effective theory, i.e. on a 
and V, and on the number of space-time dimensions. In D dimensions, using the MS 
renormalization scheme, we define 



oo 1 



s=l 

From the scale independence of the bare potentials 

v^V B = 0, (17) 
dv 

one obtains the RG equations of the different renormalized potentials. They can schemati- 
cally be written as one (vector-like) equation including all potentials: 

"■fv = B(V) , (18) 
dv 

B(V) ee - [v^Sv\ . (19) 

Note that Eq. (fl8l) implies that all the 1/e poles disappear once the derivative with respect 
to the renormalization scale is performed. This imposes some constraints on 5V: 

0(l/e) : B(V) = -2a^Z« , (20) 

0(l/e 2 ) : B (V)^ Z v + aP(a)jt.Z<p + ^^zf = , (21) 

and so on. 



6 



Out of this theory we can obtain some observables. In this paper we focus on the energy 
of the heavy quarkonium (also in its static limit). The singlet propagator near on-shell can 
be approximated to the following expression 



dte iEt d 3 R(vac\S(t,r,K)S\0,r',0)\vac) ~ n (r)0 n (r')H 



E-hf -^ B (E)+ie 



n) 



'Til 1 JYn\ 



E - El ot - 5E™ + ie 



(22) 



where n generically denotes the quantum number of the bound state: n — > (n (principal 
quantum number), I (orbital angular momentum), s (total spin), j (total angular momen- 
tum)). EP ot and (j) n {v) are the eigenvalue and eigenfunction respectively of the equation 



E p n ot U 



(23) 



and, in general, will depend on the renormalization scheme the ultrasoft computation has 
been performed with. The self-energy Y>b{E) accounts for the effects due to the ultrasoft 
scale and can be expressed in a compact form at NLO in the multipole expansion (but exact 
to any order in a) through the chromoelectric correlator. It reads (in the Euclidean) 



Vb(E) = Vi 



(D - 1)N, 



c JO 



dtre -t(K -*O r ( vac |0 E a ( t )0g.( tj 0)gE b E (0)\vac) . (24) 



The pNRQCD one-loop computation yields [231-125] 
S B (l-loop) = -^/^ 2 (l + e 



71 



2+e 



(25) 



The two-loop bare expression can be trivially deduced from the results obtained in Refs. 26- 
for the static case. It reads 



S B (2 - loop) = g 4 C f C A V*r(-3 - 4c) W(e) - (1 + 2e)vf\e) r {h B a - Ef^v , (26) 



where 



pW(e) 



(2tt) 2 4vr 2 + 2e 

1 1 
(2vr) 2 4vr 2 + 2 ' 



r 2 (l + eMe) 



r 2 (l + e) 9l (e) 



and 



2e 3 + 6e 2 + 8e + 3 2 e r(-2e - 2)T(-2e - 1) 



e (2e 2 + 5e + 3) 



(2e + 3)r(-4e-3) ' 



(27) 
(28) 

(29) 
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= 6e 3 + 17e 2 + 186 + 6 4(e + l)n / 7> 2 (e 2 + e + 1) T(-2e - 2)r(-2e - 1) 

e 2 (2e 2 + 5e + 3) e(2e + 3)N C e(2e + 3)r(-4e - 3) ' 1 ' 

From Y,b{E) it is possible to obtain 5E^ S . This will be discussed in the following sections. 

In principle we should also__ consider possible soft and ultrasoft corrections to Va- Those 
have 
Ref. 



jeen studied in Ref. |22J with LL accuracy, in Ref. [27| with NLO accuracy, and in 
7] with NLL accuracy, reaching to the conclusion that they do not contribute to the 



precision of our computation (so we can set Va = 1): 

v—V A = + 0(a 3 ) , 
dp 

whereas for the initial matching condition Va = 1 + 0(a 2 ). 



(31) 



III. V, 



A. Ultrasoft running 

We now discuss how to obtain the ultrasoft RG running of V s from T,g(E). The bare 
ultrasoft self-energy is a function of h Q — E, a and e, where this a refers to the one associated 
to the ultrasoft running. The a associated to the soft running is encoded in the matching 
coefficients (i.e. the potentials). From the combined one and two loop ultrasoft computation 
we obtain (after changing the bare alpha by the renormalized one) 



-±C f VXh -E)hh a2 M 



C f V 2 A v{h - Efv 



3'^(4tt) 2 
a{v) o, 2 (u) 



+C f V 2 r(K ~ Ef 



1087T 2 



18/3o ln z 



3^ 
_a{y)_ 

9tt 
h n -E 



367T 2 

'h, 



47 

y 



2tt 2 



10 



T F nj) 



6 In 



E 



+ 61n2-5 



6(Ca(13 + 4tt 2 ) +2/3 (5-31n2))ln 



h n -E 



-2C A (-84 + 391n2 + 47r 2 (-2 + 31n2) + 72C(3)) 



+/3 (67 + 3tt 2 - 60 In 2 + 18 In 2)) 



r + 0(e,a 3 ) 



(32) 



Note that the 1/e 2 term comes from two sources: the two loop bare result and the 1/e 
inside in the one loop ultrasoft. Quite remarkably, there is no log dependence of the 
1/e and 1/e 2 terms. This is of fundamental importance for renormalizability and a check of 
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consistency. This result has to be reexpressed in terms of the potentials of the singlet/octet 
Hamiltonian and h s — E. Positive powers of h s — E do not contribute to the energy, as they 
cancel powers of l/(h s — E) in the Green function. Therefore, we are rather interested in 
the identity (valid in D dimensions) 

,(». - Efr = W - ^ [p, [p, V?>]] +^|{P 2 , A^} + (4vr) 

1 



2m. 

+0({h s -E)) 



(AV) 2 (3d - 5) + AAV r—AV + AV 



dr 



r^-AV 

dr 



AV 



(33) 



where we have approximated h Q — h s = V} ^ — V} u> , which is enough for our precision, and 
defined AV = V o {0) - V s {0) . We have used the combination ((r-^AV) + AV), since it has a 
0(e) suppression with respect to AV. Note that in D dimensions the static potential has the 
following expansion in terms of the bare coupling constant (C s = —Cp and C Q = 1/(2N C )): 



(o) 



i/(o) _ r , a 2 9b n \ ' 

v s/o,b ~ ^s/o9b2_^ r 

n=0 



(34) 



We can now obtain the counterterms of the singlet Hamiltonian due to the ultrasoft 
divergences up to NLO in the following compact expression 



i 



2m r 
1 



(AV) 2 (3d-5) +4AV 



r—AV I + AV ) + 

dr 



r^-AV 

dr 



CfV 2 A 



a(v) a 2 (v) 



(Ca(- 



47 



2ti z 



10 



T F n } ) 



AV 



+ #p s vi-fo (47r)2 



(35) 



3tt 36tt 2 v ~ v 3 3 

The whole 1/e 2 term fulfills Eq. (12 ip and hence it is a check of the two-loop computation. 
In the above expressions we choose to keep the full ^-dependence (also in the potentials, see 
Eq. (I34p ). This is not particularly relevant for the computations we perform in this paper, 
but will be potentially important once potential divergences are included, necessary for a 
complete NNNLL evaluation of the heavy quarkonium mass. As we have already mentioned, 
in this prescription we also take vj ^ and V^ in d dimensions. Since we only need them at 
one loop, the ultrasoft divergences of Vjy°j do not show up yet. This actually means that the 
renormalized potentials are equal to their bare expressions (which one can find in Ref. 
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vj?^ ~ K/oB' anc ^ K/oB are fi n ite when we take the e — >■ limit. Moreover, at one loop, we 
also have the equality Vo°^ = —1/(N 2 — l)Vs . We prefer the scheme defined by Eq. (155]) . 
since it allows to keep the ultrasoft counterterms in a very compact manner. One is always 
free to change to a more standard MS scheme. Note that the MS scheme in momentum and 
position space are different. 

Using Eq. ( 13 5 j) and Eq. ( 12 Up we obtain the following RG equation that resums the 
ultrasoft logs of the potential: 

B Va , (36) 



where 



c s v 2 A 

1 

+ 



v\AVf + ±(aV (r^<°>) + (AV) 2 ) - -L [p, [ P ,V™]] (37) 
2a(u) a 2 (v) 



2m 2 r 



{P 2 ,A^} 



-2tx 1 



T F n f ) + 0(a^ 



3tt 9tt 2 3 " ' ' 3 

and now one could take the four- dimensional expression for the potentials. This result holds 
true in both schemes, the MS and MS (in a way this is due to the fact that the sub divergencies 
associated to a also change to make the result scheme independent). After solving the RG 
equation we find 

2 



$V S R G (r;v s ,is) 



r 2 (AVy + 



m T 



AV { r Tr V ° 0) ) +iAV)2 )) F 



1 [p, [vM Q \r)F{v s -v)]]+^- 2 { V \AV{r)F{v s -v)} 



2m 2 . 



, (38) 



where we define 



F(u s ; v) = C f V 2 



2tt 



2 h,°M 



Pa [3tt a(y„) 

Pi 1 



(39) 



3/9 (47r) 



~{a{v) -a(u,)) 
-C F Vll-\n- + 0{a 2 ) 



27tt 



(C A (47 + 6tt 2 ) - 10T F n f ) 



3 7T U a 

Note that, formally, we can set v s = 1/r as far as F(l/r,v) is kept inside the 
(anti) commutators, i.e. in the way displayed in Eq. (138]) . Yet one should be careful, as 
objects that are ill defined (not even distributions) could appear. We discuss this issue 
further below. 
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From Eq. (1381) we can easily identify the RG contribution to each potential 2 (now we 
work in the equal mass case). For the static potential we have 

SV s %(r; u„ u) = r 2 (AV) 3 F(u s ; p) , (40) 

which agrees with Ref. Q. 

For the 1/m potential the RG contribution reads 

W s %(r; u„ u) = a(aV (r±V^ + (AV) 2 )f(i/.; p) . (41) 

For the momentum- dependent 1/m 2 potential we have 

SV ^rg(^ "s, v) = 4AV{r)F{v a - v) . (42) 

The RG correction to V r deserves a special discussion. From Eq. (!38|) we have 

W r %(r; p., v) = -2 [p, [p, V^F(v s ; u)]] = J ^e^5V r %( g] v., v) , (43) 

where 

W r %(q; p s , p) = -2 q 2 V}°\q)F(u a ; u) , (44) 

and Vo^ is the Fourier transform of Vo°\ Note that we have chosen to write 5V r in mo- 
mentum space as well. This we could also do for the previous potentials but now will be 
particularly convenient. The reason is that for SV r! RG{ r ', v a , v) we cannot naively set v s = 1/r 
(as we did for the previous potentials and was convenient to resum logs of r) because diver- 
gent distributions like 5^ 3 ^(r)lnr could appear. On the other hand we would like to make 
something analogous, as it allows to resum (and keep in a compact way) some subleading 
logs. We then choose to set p s = q in the Fourier transform and define 

SV%(r; v) = J ^e^5V%(q; q, v) , (45) 

where SV^ G (q; q, p) is Taylor expanded in powers of In q. We will elaborate on this expression 
in the following subsection (in practice we will only need the single log of this expansion). 
Finally also note that SV^ G vanishes in the large iV c limit. 

2 And from previous equations one could also easily obtain the counterterms and RG equations for each 
potential. 
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B. Initial matching conditions 



The only thing left to obtain is the initial matching condition for the potential. We 
will only consider the spin-independent terms and stick to the equal mass case, since the 
matching coefficients at one loop are not known in the non-equal mass case, although many 
partial results exist. 

We now study each potential separately. Note that the result will depend on the basis 
of potentials used and on field redefinitions except for the singlet static potential, which is 
unambiguous 3 . 

The initial matching conditions for the static potential in our MS scheme reads 

V&r;,) = -e^M{l + t(^)"^r)}, (46) 

n=l V ' 

with coefficients 

ai(v,r) = ai + 2(3 In [ye lB r) , 

a 2 (u, r) = a 2 + y /3 2 + ( Aa,^ + 2ft) In {ve^r) + 4/3 2 In 2 (ve™r) , 

571-2 

a 3 {v, r) = a 3 + a^ 2 + — ftft + 16Cs/3 3 



+ (W 2 /3 3 + 6a 2 /3 + 4 ai ft + 2/3 2 + yC> 2 ) In ' 
+ (l2ai/3 2 + 10/5 /?i) In 2 (z/e 7B r) + 8/3 3 In 3 (z/e 7E r) . (47) 

Explicit expression for a, can be found in the literature [l, 2, 22, 24, 25, 29, 30]. For ease of 
reference we display them in the Appendix. 



[ve lE r) 



We choose to write the initial matching conditions of the potentials, like Eq. f f46|) . in terms 
of the (single) factorization scale of the effective theory v and not u s , as u s does not appear 
in a matching computation done order by order in a. The (left-over) factorization scale 
dependence of the matching coefficient would cancel with the scale dependence of loops in 
the effective theory (in practice most of the scale dependence in the potentials, for instance 
in Eq. (115]). cancels with the scale dependence of a(v), effectively becoming a(l/r) but 
this does not change the physical principle). When a RG analysis is done one introduces 



This observation also applies to the RG corrections to the potentials obtained in the previous subsection. 



12 



two factorization scales: the one where the running starts, which we have named v s in the 
previous section, and the factorization scale that would cancel with the scale dependence 
of loops in the effective theory, which we have always named v. Therefore, the total RG 
improved static potential then reads 

Kw™^ v ) = V lU r > ^ + 6V ^c(r; v s , v) , (48) 

and it is correct with NNNLL accuracy. 

For the 1/m potential the initial matching condition reads 

'bxPa 2 



s,ms v ' > 2r 2 
where 















, b 2 = - 


89 2 
~36 ^ 



2 



- 2 -(Cl + 2C A C f )"j In^Ve 2 ^ 



(49) 



Ct 89 1 7 49 

h = -£-C A , b 2 = --C\ + —C A C f + -C A T F n f - ^C f T F n f . (50) 



The two loop result has been taken from Ref. [3jJ and changed accordingly to fit our 
renormalization scheme for the ultrasoft computation. This explains the different coefficient 
6 2 we have compared with that reference. Ours is the proper one to be combined with the 
ultrasoft correction obtained in Eq. ( \72\i in the next section. The total RG improved 1/m 
potential then reads 

v lm G ^ v ) = v !U r > + 6V sUr; v s , v) , (51) 

and it is correct with NNLL accuracy. 

For the momentum-dependent 1/m 2 potential the matching coefficient reads at one loop 



(52) 

The total RG improved momentum-dependent 1/m 2 potential then reads 

v Sm^ ") = v Sm^ ^ + 5V $™(r; ») , (53) 

and it is correct with NLL accuracy. 

The results obtained for the above potentials are exact to the required accuracy, since 
their soft running is trivial and there are not potential loops at this order. This can be 
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traced back to the fact that there is no dependence on the matching coefficients inherited 
from NRQCD, and no logarithms proportional to the mass appear. The independence of the 
potential on u s , v s ^r — 0, can be easily implemented by setting u s = 1/r, giving information 
on the In r dependence. In fixed-order bound state computations is convenient to work using 
a(v s ) as the expansion parameter with v s ~ ma. By expanding a(l/r) in powers of a{y s ) 
times ln(z/ s r) subleading corrections appear that have to included with the appropriated 
precision. 

The case of Vr 2 ^ is different. It depends on the mass of the heavy quark (i.e. on the 
matching coefficients inherited from NRQCD). In this case we can not give the complete 
expression with NLL accuracy. For such accuracy we would need the soft divergences to 
a higher order and the inclusion of effects due to potential divergences. The latter are 
absorbed in delta-type potentials. Note that when computing their running at NLL the 
ultrasoft LL running would enter indirectly. Such kind of computations have already been 
undertaken in Refs. (5), |(3] for the spin-dependent corrections to the hyperfine splitting of 
heavy quarkonium, a similar analysis for the spin-independent corrections would go much 
beyond the aim of this work. Either way we can still give the NLO matching coefficient for 
the potential, the relevant starting point for a complete NLL. It is first convenient to write 
it in terms of the potential in momentum space 

d 3 q 



r,MS V ' > 



(27T)" 



(54) 



where 



r,MS 



a{q){l + c D {v)-2c 2 F {v)) 



(55) 



1 1 ~ 

H — (d vs (u) + 3d vv (v) + —(d as (u) + 3d vv (u))) + 5V soft (v, q) 

71 c f 



soft 



Ol 
7Y 



9 25, v 
- H In - 

4 6 q : 



, 1 7 , v 

( 1 ^-3^7 



(56) 



The soft one-loop result 5V so f t is taken from Ref. 3]. The rest corresponds to the hard 
contribution. We have checked that upon expanding the NRQCD matching coefficients in 
powers of a(v s ) we agree with that reference at 0(a 2 ). In order to use SV so f t we hac 



change the scheme used for the Pau 
(see also the discussion in Refs. 



i cr matrices in the computation of d vv in Ref. 



to 
3% 



33]). The new expression for d vv can be found in the 
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Appendix, as well as for the other NRQCD matching coefficients. Eq. is correct at LO, 
LL (soft). It is also correct at NLO, providing with the right initial matching condition to 
compute the NLL result. We can then give the following expression for the potential 



V 



(2),RG/ 



(r; v) = V^(r; v.) + 5V® {r; ^ „) . 



(57) 



r,MS v ' > r.MS 1 

This expression does not incorporate the NLL running associated to the Fourier transform 
of lng yet. In order to do so it is rather more convenient to consider 



where 



r,MS 



V^ RG (r; u) = V r %(r) + 5V}%( r] u) 



d 3 q 



3 iq-ri 



V {2 -L(a a) 



(5? 



(59) 



'> MS J (27r) 3 r ' MS 

Although this new expression, Eq. ( 1581) . does not account for all the logs at NLL order, yet 
it does for the ultrasoft logs and those associated to In q. This is enough for the computation 
of the mass of the / ^ and s = states with NNNLL accuracy. The reason is that the 
delta potential does not contribute to the mass of those states. This is more clearly seen if 
we Taylor expand the above expression in powers of (the Fourier transform of) In q: 



V 



(2),RG 



r,MS 



(r; v) 



d ,~ 



(In v s )qj-(V^(q, q) + 5V r %{ q] q, u))\ q ^ 



(2) 



dq 



where we have used 



47T \ r 6 



d 3 q 



In q . 



(60) 



(61) 



Only the term proportional to reg^ contributes to the mass of the states with I ^ 0. Higher 
order terms in the Taylor expansion are subleading. The derivative of the potential can be 
basically written in terms of the NRQCD matching coefficients 



dq 



C s a\v s ) 



^T F n f (c D + c?) + {f3 -^C A )c 
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- y + (yQ - -C A )c k 



(62) 



q-^( S V}$ G {q;q,v))\ q=Us 



1 + ln 



(63) 
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In order to obtain Eq. (ISTJI) . we have used the explicit a(r _1 ) dependence of the potential. 
At the end we have a double expansion in a(v s ) and (x(y). In the bound state calculation 
they get replaced by v s ~ ma and v ~ ma 2 . 

We take Eq. ( 1581) or Eq. ( 1601) as our final expressions for V ^ (r; z/). 

We have then obtained the RG improved and initial matching coefficients for the different 
potentials that compose V s to the order of interest. As we have already mentioned the 1/m 
and 1/m 2 potentials suffer from field redefinitions ambiguities. Therefore, in some circum- 
stances it can be convenient to cast the initial matching conditions of V s in the following 
unified form 4 : 

W; «.) = Cfr; ,.) + + ^(i{ P ». vSW" «"> } + • («) 

as well as the total potential to be introduced in the Schroedinger equation 
V s%^ v ) = V sM r 'i v * = + SV s>RG (r; vs = 1/r, i/) (65) 

y(l),flG/ \ 

Note that we have worked in a basis where Vjji is set to zero, at least within the accuracy of 
our computation. 

IV. OBSERVABLES 

A. Static potential and Energy 

The singlet static energy can be considered to be an observable for our purposes. It 
does not suffer from field redefinitions ambiguities and can be easily compared with the 
lattice determination. It consists of the static potential, which is a Wilson coefficient, and 
an ultrasoft contribution, both of them can be taken either bare or renormalized. More 
explicitly, the expression for the static energy with NNNLL precision reads 

E ^ = V lm G ^ » = W) + 8E®£{r- t v = AV) , (66) 



4 Note that we use V\^(r) and not V SLfr, v s 

rMS y ' r,MS y 
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where 



<W M) = -Cf^Vfvf^- ( 6 ln^ + 61n2-5) . (67) 

Note that one has to be careful in using the same scheme in the soft and ultrasoft to obtain 
the correct result for the energy of the static singlet state. For ease of reference we also give 
the NNNLO result. It reads 

E NNNLO {r) = y(V_ {r . y) + SE ^(r; v) . (68) 



We can also give the expression for the subleading ultrasoft contribution to the static energy: 
<^Vo = C f AAVfV^ (W„ln= (69) 
-6 (C A (13 + 4tt 2 ) - 2/3 (-5 + 3 In 2)) In (^p) 

-2C A (-84 + 39 In 2 + 4vr 2 (-2 + 3 In 2) + 72((3)) + f3 (67 + 3tt 2 - 60 In 2 + 18 In 2 2) J . 

This contribution is relevant for the complete N 4 LO or N 4 LL computation of the static 
energy. For the former the only computation left is the soft four-loop contribution to the 
static potential. 



B. Energy l/0s = 

We are now in the position to obtain the heavy quarkonium energy with N 3 LL accuracy 
(for I j£ and s = 0): 

E n l js \ ¥0 ,s=0 = ( E ^ nljs + 5E WS,nl) li#0,s=0 , (70) 

where , . is the eigenvalue of the equation 

MS,nljs to H 

{jk + V Ss) = E mn^s(r) , (71) 

and V s % is Eq. ([65]) minus + ^) p 4 . 

The exact solution of Eq. (1TT]) correctly produces all necessary soft and potential terms 
for the aimed N 3 LL accuracy, as well as some subleading terms. Such exact solution would 
only be possible to obtain through numerical methods (which on the other hand could 
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actually be more easy to implement in practice). If we want to restrict ourselves to a strict 
N 3 LL computation (in particular if seeking for an explicit analytical result), Eq. flTTj) should 
be computed within quantum mechanics perturbation theory up to NNNLO for general 
quantum numbers. Up to NNLO such computation was performed in Ref. 34j. The lengthy 
N 3 LO computation is missing, beyond the aim of this work, and will be considered elsewhere 
(for I = and n = 1 such computation has been performed in Ref. 4J). 

The ultrasoft correction to the energy due to the ultrasoft correction can be written in 
the following compact form 



U£/ MS,n.i 



a 



(n,l\[C f r(h -E n>l ) 3 V2 



(f>ln p^)+61n2-5 



(72) 



1087T 2 



1800 In 2 



h„ - E, 



n.l 



- 6 (N c (13 + 4tt 2 ) -2/3 (-5 + 31n2))ln 



h Q — E n i 



v 



+2C A (84- 391n2 + 47T 2 (2 -3 In 2) - 72((3)) + (3 (67 + 3vr 2 - 60 In 2 + 18 In 2 2))] rj\n,l) 
where the states |n, I) and the energies E n j used above are the solution of the Schroedinger 



potential including the 1-loop static potential (i.e. with NLO accuracy): 
^ 2 C f a s (u) 



P 



2m r r I Ait 



ax(u; r) 



\n,l) = E n} i\n,l) . 



h could also be approximated to its NLO expression: 



h a 



1 at.Jv) 



l + -^ai(i/;r)k 



(73) 



(74) 



2m r 2N C r I 4tt 
Eq. ( I72|) includes the complete LO (D(ma 5 ) and NLO 0(ma 6 ) ultrasoft effects, as well 
as subleading effects. The LO expression would be enough for the N 3 LL precision. A 



semianalytic expression exists for / = states [25] but missing for general quantum numbers 
and would require a dedicated study, again beyond the aim of this paper. In a strict fixed- 
order computation one should expand the wave functions to the appropriate order, as well 
as h — E n j, but in some situations it could be more convenient to handle this expression 
numerically. 



V. CONCLUSIONS 



We have computed the NLL ultrasoft running of the spin-independent singlet potentials 
up to 0(1/ m 2 ), and the corresponding contribution to the spectrum. This includes the static 
energy at NNNLL order. 
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Whereas there are still some pieces left for a complete NNNLL evaluation of the heavy 
quarkonium spectrum, our result provides with the missing link for the complete result for 
I 7^ and s = states, the structure of which is given in this paper for the first time. The 
full explicit analytic form will be presented elsewhere. 

Another by-product of our computation is the N 4 LO ultrasoft spin-independent contri- 
bution to the heavy quarkonium mass and static energy. 

Note added: After our paper appeared on the web, the preprint 



35| was sent to the 

arXives. In this reference expressions for the NLL ultrasoft running of the 1/m and 1/m 2 
potentials were given in the vNRQCD framework and agreement with our results claimed. 
Note that in order the expressions for the 1/m 2 potentials to be equal, the contribution 
associated to Eq. 31 in 35] has to be added to the result obtained in Ref. jl4j |. It also 



remains to be explained how (and if) the dependence on the infrared regulator mentioned 
in Ref. has disappeared. 
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Appendix A: Constants 



7>4 C A = N C - C f =^. (Al) 



0o = H% - \T F n f - (3, = 34^ - ^-C A T F n f - AC f T F n f - (A2) 



A = ^C\- l ^C^ A T F n f ^C A nn}-^ (A3) 

31CU-20TP7V .... 

«i = g -; (A4) 



400 V T F 2 /55 \ 

a 2 = ^ C f rif T F I— - 16 C(3) J 



/4343 16 7r 2 -vr 4 22((3)\ /1798 56 C(3) 



+ 
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n - „(3)„3 i „( 2 )„2 , (1) , (0) 



(A5) 



where 



,(3) 



(2) 



/ 12541 368C(3) 64tt 



243 + 3 
709.717) C\T F 



+ 



/ 14002 416C(3) 



135 J V Q1 :? 

71281 



162 



+ 264C(3) + 80C(5) C A C>T F 



(;'« d^abcd 



/286 296C(3) . A ^ 9 , . . . 

( — + — ^ - 160C(5)J C)T F + (-56.83(1)) 



,(°) 



502.24(1) C\ - 136.39(12) 



^abcd^abcd 

N A 



and 



d abcd d abcd N C (N 2 C + 6) 



(A6) 



(A7) 



A^4 48 

The NRQCD matching coefficients have been computed over the years. In Ref. [36| one 
can find the NRQCD matching coefficients of the one heavy quark sector at NLO, whereas 
in Ref. 32| one can find them for the two heavy quark sector (although with a difference of 



scheme for dm,). In Ref. [37j the LL soft running of the one heavy quark sector was obtained 
and in Ref. [38[ the LL soft running for the two heavy quark sector. a(m) has n/ active 
light flavours: 



d 2 (m) ■■ 
c F (m) = 1 + 
co(ra) = 1 + 



a (to) 
" 60tt 
a(rn 



-Tv 



2ty 
a(rn 



2tt 



(C f + C A ), 
Ca — 16d 2 (w). 



(A8) 
(A9) 
(A10) 



d a ss (m, 
d a 



sv 



m 



d a vs (m) 



d a vv {m) 



C f ) (2-21n2 + m) 



a {m)Cf 
0, 

^(- 3 ~C A + 4C f )(2-2ln2 + in) 



-mm 



aim, , 
1 + ( T, 



7T 



R 



1 ^ , 5 

— rif 2 hi 2 Z7r 

3 S \ 3 



(All) 
(A12) 
(A13) 



36 ' 



(A14) 
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d ss (m) 
d sv (m) 
d vs (m) 
d vv (m) 



d a ss (m) 3d a sv (m) N 2 - 1 N 2 - 1 2_ 



— C/ I a 2 (m) , 



<£„M N 2 C - 1 ja . . iV c 2 - 1 



2iV„ 



+ 



2iV„ 



47V 2 vs 

c 



— C/ 1 a (m) , 



-C(^) - 3^(m) + + 3 ^ (m: 



2/X, 



2N r 



+ + a 2 (m) , 



-C(-) + <fc(m) + ^ - *M + a 2 (m)(2C f - ^) . 



(A15) 



We now define z 



a(m) 



00 



l/(27r)a(i/.)lii(£), 



c F (u s ) = c F (m) - 1 + z~ Ca , 
cd{v s ) = c D (m) - 1 + 



d ss {v s ) 

dsv(Vs) 

d vs {v s ) 



9C A f 5C A + 4T F n f _ 2Ca C a + 16C f - 8T F n f 
9C A + 8T F n f \4C A + 4T F n f Z + 2(C A - 2T F n f ) 



-7C A + 32C A C f - 4C A T F n f + 32C f T F n f 4Tpnt /3 _ 2 c A /3 
A(C A + T F n f )(2T F n f -C A ) 



8T F rif 
9C A 



-2C A , (20 32 C> \ 
\13 13 C A J 



d ss (m) + 4C, (c, - ^fj ja{m) [z*> - l] , 
d sv {m) , 



8vr 



d vs (m) - (C f - C A ) —a{m) [z? - l] 
Po 



27C 2 A 



IT 



9C A + 8T F n f p t 



-mm) 



5C A + 4T F n f Po 
~4C A + 4T F n f f3 -2C A 



Po-2C A 



1) 



C A + 16C f -8T F n f „ 

2(C A -2T F n f ) { > 
-7C 2 A + 32CaC/ - 4C A T F n f + 32C f T F n f 
+ 4{C A + T F n f ){2T F n f -C A ) 



x- 



3/?o 



3p + AT F n f - 2C A 



y /3 +4T F n f /3-2C A /3 



i) 



ITpUf 

~9C A ~ 

x 



Po 



d vv {v s ) = d vv (m) + 



Po - 2C A 
[zP° - 1] - 



■>*( 



6/3 



6/3 - 13C A 



20 32 C/ 
13 + 13 

- 1 



-am) {z^- 2Ca - 1} . 



(A16) 



\ — 2C A 

The results displayed above for the NRQCD matching coefficients, c's and ofs, are correct 
with LL and NLO accuracy, but not beyond. 
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We will also need 



mi \ = 9C A f 5CU + 47>n/ 2C/ 
1 1 s> ' ' 9C A + STpUf \ AC A + ATpUf 



C A + 16C/ - 82>n/ 
2(CU-27>n / ) 



-7C| + 32C A C f - AC A T F n f + 32C f T F n f 4 r F n f /3-2C A /3 
4(C A + T F n / )(2T F n / -C A ) 



20 
13 



32 Cj_ 
13 



1 - z" 



(A17) 



This is correct with LL accuracy, which is enough for us. 
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